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Resumo. Neste trabalho fizemos comparagoes entre formulagoes semi-discretas 
para a obtengao de solugoes numericas para a equagao ID de Burgers. As formu¬ 
lagoes conslstem em discretlzar o domlnio temporal via metodos impllcitos multi- 
estagios de segunda e quarta ordem: aproximantes de Fade Rn e R22', e o domlnio 
espacial via metodos de elementos finitos: mlnimos quadrados (MEFMQ), Galerkin 
(MEFG) e Streamline-Upwind Petrov-Galerkin (SUPG). Conhecendo as solugoes 
anallticas da equagao ID de Burgues, para diferentes condigoes iniciais e de fron- 
teira, foram realizadas analises dos erros numerlcos a partir das normas L2 e L^a- 
Verificamos que o metodo com o aproximante de Fade R22 adicionado as formu¬ 
lagoes MEFMQ, MEFG e SUPG, aumentou a reglao de convergencia das solugoes 
numericas e apresentou maior precisao quando comparado as solugoes obtidas por 
melo do aproximante de Fade Rn. Gonstatamos que o metodo R22 amenizou as 
oscilagoes das solugoes numericas associadas as formulagoes MEFG e SUPG. 

Palavras-chave. Equagao de Burgers, Aproximantes de Fade, Metodos impllcitos 
multi-estagios, Metodos de elementos finitos. 


1. Introdugao 

Com a evolugao da mecanica computacional ocorreu a intensificagao de pesquisas na 
resolugao numerica de equagoes diferenciais com impactos positivos para a sociedade 
m- Devido a vasta aplicabilidade destas pesquisas em problemas que envolvem 
processes convectivos e difusivos miniiiniiiiiisiiii], interessa-nos as formulagoes 
numericas que possam ser aplicadas a estes problemas, em especial a equagao de 
Burgers O ^ [3 [HI [13 US [21 US] . 

Varies autores apresentaram solugoes numericas para a equagao ID de Burgers, 
usando metodos de elementos finitos |3[Z1[H1[IS1[I1[IS1[2III21I2S], de elementos 
de contorno de diferengas finitas IS [H [H], assim como alternativas para 
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resolver o termo nao linear convectivo como esquemas upwind jO] [10] e tecnicas 
de linearizagao muini. Neste contexto faremos comparagoes de formulagoes 
semi-discretas para a equagao ID de Burgers, onde utilizaremos formulagoes semi- 
discretas para a discretizagao temporal e espacial e realizaremos uma linearizagao 
no termo convectivo. Esta linearizagao altera o tamanho do elemento em cada etapa 
usando a informagao a partir do passo anterior IS 1131116] transformando a equagao 
de Burgers em um problema linear local. 

As formulagoes semi-discretas consistem em discretizar o dominio temporal uti- 
lizando metodos imph'citos multi-estagios e o dominio espacial via metodos de ele- 
mentos finitos [HllMlIlS]. Em particular, para resolver a equagao ID de Burgers, con- 
sideramos os metodos imph'citos multi-estagios de segunda ordem, Rn, e de quarta 
ordem, R 22 , e tres formulagoes de metodos de elementos finitos: minimos quadrados 
(MEFMQ), Galerkin (MEFG) e Streamline-Upwind Petrov-Galerkin (SUPG). 

Gomo o objetivo deste trabalho e, conhecido a solugao analitica da equagao 
ID de Burgers, determinar quais formulagoes semi-discretas apresentam solugoes 
numericas com precisao em um intervalo maior de tempo, aumentando desta forma 
a regiao de convergencia, apresentamos resultados comparatives entre as solugoes e 
simulagoes numericas para diferentes tempos. Apresentamos tambem tabelas com 
os calculos dos logaritmos dos erros entre as formulagoes utilizando as normas L 2 e 
Loo • 

Desta forma, este artigo encontra-se estruturado como segue: a segao 2 apresenta 
a equagao ID de Burgues na forma geral e tambem na forma padrao de operador; 
na segao 3 e apresentada a discretizagao espacial obtida por meio dos aproximantes 
de Pade; na segao 4 sao apresentadas as formulagoes de elementos finitos, utilizadas 
na discretizagao espacial da equagao ID de Burgues; na segao 5 e introduzida a 
tecnica de linearizagao do termo convectivo transformando a equagao de Burgers 
em um problema linear local; a segao 6 mostra os resultados numericos obtidos no 
trabalho e finalmente as conclusoes sao apresentadas na segao 7. 


2. Equagao ID de Burgers 


Seja a equagao ID de Burgers definida em um dominio = [a, b] G R limitado e 
aberto com fronteira P = 9D, satisfazendo 


Ut{x,t)u{x,t)ux{x,t) - eUxx{x,t) = f{x,t) em D, (2.1) 

u{x,0)=Uq(x) y X G it, (2.2) 

u{a,t) = Co = u{b,t) em r, (2.3) 


onde e = 7?e e o numero adimensional de Reynolds; u{x,t) e a componente 

da velocidade do fluido na diregao do eixo x e f{x,t) representa o termo fonte. A 
equagao (2.2) e uma condigao inicial com uo sendo uma fungao dada e (2.3) e uma 
condigao de fronteira do tipo Dirichlet, sendo cq uma constante. 

Reescrevendo (2.1 ) na forma padrao de operador, temos 


ut -f C{u) = /, 


(2.4) 
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onde C{u) = uux — tUxx representa o operador espacial e descreve a soma dos 
operadores nao linear e linear, convectivo e difusivo, respectivamente. 


3. Discretizagao temporal 

Muitas tecnicas numericas para a discretizagao temporal sao utilizadas para resolver 
equagoes diferencias parciais |8l|9l|23]. Consideramos a tecnica passo de tempo 
obtida por meio dos aproximantes de Fade, cuja aproximagao para o operador de 
evolugao e dada por 


£{M) :w(r) 


(3.5) 


que permite transportar a solugao numerica em um determinado tempo = nAt 
para o proximo tempo + At, sendo At um passo de tempo. Nesta situagao, 

a evolugao de £(At), obtida a partir do desenvolvimento da serie de Taylor, e dada 
por meio do operador exponencial, on seja 


u(t"+i) 


1 .9 


d 


1 9” 

... + -Ar— 
n\ df^ 


exp ( At^ ) u{t^). 


(3.6) 


Assim, a tecnica de passo de tempo de varias ordens pode ser obtida usando os 
aproximantes de Fade para o operador exponencial e^, sendo h — At^ [5] [S]. Os 
primeiros aproximantes de Fade para encontram-se apresentados na Tabela 1. 


Tabela 1: Aproximantes de Fade da fungao exponencial 


Rlm 

M=0 

M=1 

M=2 

M=3 

L=0 

1 

1 

‘A 

b 

1-h 

2-2h+h‘^ 

6-6h+3h‘^-h^ 

L=l 

r+7t 

2+Ti 

Wp2R 

:24+6/i 

1 

2-h 

6-4h+h^ 

24-18h+6h^-h^ 

L=2 

2+2h+A 

6+4h+A 

12+6h+fd 

60+24)1+37? 

2 

6-2/1 

12-6h+h^ 

60-36/i-9/i2-h3 

L=3 

6+6h+3h^+h^ 

6 

24+lSh+6h^ + h^ 

24-6/1 

60+36h+9h^+h^ 

60-24/i+3?i2 

120+60/1+12/1^+/!^’ 

120-60/i+12/i2-/i3 


Em HU EH] os autores mostram que a aproximagao de Fade Rlm e incondici- 
onalmente estavel se satisfaz a condigao M — 2 < L < M. Assim, as tecnicas de 
multi-estagios imph'citas empregadas na parte da discretizagao temporal, utilizadas 
neste trabalho, sao incondicionalmente estaveis com erro de truncamento de ordem 
C>(At2") [5]. Logo, podemos reescrever o metodo implicito na seguinte forma de 
operador 


— - WAut = w<. 

At * 


(3.7) 
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onde o vetor Au tern n componentes (ou 'figsiagios) 0 ® ® ^ derivada parcial 

de Am com respeito ao tempo. Substituindo \2A\ em e considerando que o 
operador £ e linear com coeficientes constantes, temos que 

/\ 7 / 

—+W£(AM)=w[r-£(M")]+WA/ (3.8) 

com Am, W, A/, e w dependentes do metodo escolhido. Para os metodos impHcitos 
i?ii e i? 22 ) os valores de Am, W, A/, e w encontram-se definidos nas segoes 3.1 
e 3.2. Devido a limitagao de espago, definimos apenas as formulagoes compactas 
destes metodos, outros detalhes podem ser obtidos em iBtni- 


3.1. Metodo de segunda ordem 

A formulagao compacta i?ii, Crank-Nicolson 0E], e dada por 

Am = m"+i - m", a/ = /"+! - /" 
W = l/2, w = l. 


3.2. Metodo de quarta ordem 

A formulagao compacta R 22 [Hlin] e 


r m”+i/2 - m" 1 
I m”+i - m"+i/2 J ’ 



7 -1 

13 5 ’ 


A/ = 


^* 71 + 1/2 _ jn 



(3.9) 


(3.10) 


4. Discretizagao espacial 

Neste trabalho utilizamos o metodo de elementos finitos para a discretizagao es¬ 
pacial. Este metodo introduz fungoes bases {^po ,..., g?m} que geram o subespago 
onde esta sendo procurada a solugao exata, com suporte localizado nos pontos no- 
dais dos elementos El- Para definir as fungoes bases realizamos uma discretizagao 
no intervalo [a, 6], dividindo-o em m sub-intervalos (ou elementos) Cj = (xj-i, Xj), 
j = 0,1,...,TO, de comprimento hj = Xj — Xj-i. 

Seja Uh uma fungao teste linear em cada elemento e_,- e conti'nua sobre [a, h] 
satisfazendo as condigoes de fronteira Uh{a) = cq = Uh{b). Escolhendo os valores 
Mo, Ml, M 2 , ..., Um nos nos Xj, como parametros para descrever Uh{x) e expressando- 
os sobre cada elemento Cj, obtemos Uh{x) = 'ij}i\x)uj-i + il) 2 \x)uj com x € Cj, 
onde 

V 4^^ / \ ) ' 


(x) 


(4.11) 
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Podemos escrever Uh{x), sobre todo o dommio [a, 6], como sendo 

Uh{x) = ipoix)uo + (pi{x)ui + . . . + ipmix)Um, (4-12) 

onde 


1 


X G Cj; 

ipj(x) = < 

1 ^i^^\x), 

X G Cj + i 1 < J < TO — 

\ 

[ 0, 

caso contrario; 

Foix) = <1 

'ip[^\x), 

a; G ei; 

. 0, 

caso contrario; 

iPn{x) = < 

f 

1 0, 

X G Cm , 
caso contrario; 


sao as fungoes bases e ipj (x) e uma fungao linear seccionalmente conti'nua, com valor 
um para o no Xj e nulo para os outros nos. Portanto, Uh{xj) = Uj. Especificamente, 
se temos uq = Cq entao (4.12) satisfaz as condigoes de fronteira ( |2.3| ) e para os 
parametros ui,U 2 , ■ ■ ■, Um valores arbitrarios. Todas as fungoes Uh{x) constituem o 
espago das fungoes teste 14. Considere o espago de Hilbert V = onde V e 

o conjunto das solugoes tentativas. Assim, seja 14 o subespago de dimensao finita 
do espago de dimensao infinita V, formado por fungoes lineares seccionais geradoras 
de um conjunto de m elementos de V denotado por 14 = ['-Po, ■ ■ ■, <Pm\- As fungoes 
bases ipj sao obtidas a partir do metodo de elementos finitos, considerando a partigao 
o = 2:0 < a;i < a ;2 ... < Xm-i < Xm =b. 


4.1. Metodo de elementos finitos via mimimos quadrados 

O metodo de elementos finitos via mimimos quadrados (MEFMQ) consiste em apro- 
ximar os termos da parte espacial de (13 por meio de uma formulagao variacional 
El, obtida usando o metodo de minimos quadrados que minimiza o quadrado da 
integral do residuo. 

Seja V o espago de Hilbert e defina o funcional 

F :V ^ R 

(4.13) 

onde F(u) = j^[u(x)ux(x) — tUxx(x) — f{x)Ydx para todo x G ft sobre todos os 
u G V. Minizando F com respeito a para j = 0,1; 2,..., to, e usando a 
derivada de Gateaux [1], podemos resolver o seguinte problema variacional: achar 
£ V tal que 

= Fm(u;w), Vw G P, (4.14) 

onde definimos o funcional bilinear : V xV ^ R e o funcional linear 

Fm{u] ■) :V ^ R por [T7] 

/ [e^UxxWxx - euUxWxx]dx, \/w gV; 

Jq 


aM{u-,Ux,w) 


(4.15) 



6 


Ladeia, Romeiro, Natti e Cirilo 


Fm{u;w)= / [-€w^j;f+ Uxwf+ uw^f]dx, Mw &V. (4. 

Jn 


16) 


Obtido a formulagao variacional, podemos resolver a parte espacial de (2.1) utili- 
zando MEFMQ, considerando o subespago 14 C E, que consiste em: determinar 
uma solugao aproximada G 14 tal que 


= FM{uh\Wh), ^Wh G V 


(4.17) 


onde V = Ux- 


4.2. Metodo elementos finitos via Galerkin 


Uma segunda aproximagao para a parte espacial de (2.1) e realizada atraves do 
metodo de Galerkin (MEFG) |2I1]- Seja V o espago de Hilbert e dada a formulagao 
fraca, podemos resolver o seguinte problema: determinar G V tal que 




7) = Fciw), Vw G V, 


(4.18) 


onde definimos o funcional bilinear aciu;-,-) : V x U —>■ R e o funcional linear 
-F’g(’) -V^R por |T7] : 


aG{u-,Ux,w) 


/ [uuxw - eUxWxx]dx, \/w GV; 
Jn 


(4.19) 


Fg{w) 


/ wfdx, Vre G V. 

Jn 


(4.20) 


Para resolver a parte espacial de (13, utilizando MEFG, consideramos I 4 C V, no 
qual o problema agora consiste em determinar uma solugao aproximada G 14 
tal que 


= Fciwh), ywh G 14 . (4.21) 


4.3. Metodo estabilizado Streamline- Upwind Petrov-Galerkin 

O metodo estabilizado Streamline- Upwind Petrov-Galerkin (SUPG) contorna as li- 
mitagoes do metodo de Galerkin [5] . O SUPG e uma combinagao da formulagao de 
Galerkin com termos baseados no resi'duo, em m'vel de elementos. Estes termos sao 
balanceados por parametros de estabilizagao, resultando em formulagoes variacio- 
nais consistentes com as propriedades de estabilizagao, superiores as da aproximagao 
de Galerkin El [7]. 

Assim, o metodo estabilizado SUPG para aproximar a parte espacial de (2.1), 
consiste em determinar Uh G 14 tal que 


aoiuh; Vh, Wh) + Esupoiuh; Vh, Wh) = Fciwu), ^Wh G 14, (4.22) 

onde = Vh e EgijpG{uh',Vh,Wh) indicam os termos de perturbagao adicionados 
a formulagao variacional padrao (|4.21|). Estes termos sao adicionados de forma a 
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preservar a consistencia do metodo para obter a estabilidade numerica, dada pela 
expressao 


EsuPG{uh',Vh,Wh) 



'^{V{uh] Wh), Tn{uh))ai 


(4.23) 


sendo V{w) = a perturbagao da fun^ao teste, enquanto o termo residual TZ 

e o parametro r sao definidos, respectivamente por [7]: 


d'^Uh 
^ dx"^ 


(4.24) 

4ey^ 


(4.25) 


Para resolver a parte espacial do problema (13 utilizando SUPG, consideramos 
Vh G V para n = 0,1, 2,..., iV, que consiste em determinar uma solugao aproximada 


S Vh tal que 


+ EsuPGdi^^;vl'^^,Wh) = Fdwh), 'dwh G 14- (4.26) 


5. Linearizagao do termo convectivo 


Varias tecnicas para resolver o termo convectivo da equagao de Burgers podem ser 
encontradas na literatura 0 HI 13 HSl US]. Neste trabalho, realizamos uma lineari¬ 
zagao no termo convectivo de (I13> que consiste em alterar o tamanho do elemento 
em cada etapa utilizando a informagao do passo anterior Emus], transformando 
a equagao de Burgers em um problema linear local. Para isto, multiplicamos ambos 


os lados de (2.11 por uma fungao teste w G V e integramos, resultando em 


{ut + uux — euxx — f)wdx = 0. 


(5.27) 


Como a solugao de (2.11-(2.21 e procurada sobre o dommio a < x < b, com condigoes 


de fronteira em x = a e em x = b, consideramos o subespago de dimensao finita 
Vh C V, onde as fungoes bases (pj sao obtidas a partir do metodo elementos finitos, 
utilizando a partigao de tamanho hj = Xj — Xj-i, mapeada por uma coordenada 
local cr, onde Xj = Xj-i -I- ahj, 0 < cr < 1 0. Logo, constrmmos uma fungao teste 
Uh e escolhemos os valores uq, ui, U 2 ..., Um nos nos Xj. Podemos entao reecrever a 


equagao (5.27) como 


E"./„( 


duj 

dt 


1 


V) . 

dx ‘ 


-C 


d^ipj (x 
dx^ 


-iC-v)fuj'^V>^ix)dx = 0, (5.28) 
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V ipi{x),ipj{x) S Vh, onde V = C = '^ sao localmente constantes sobre cada 

elemento e Wh = ‘fii{x), para i = 0,1, 2 ,..., to, e definida como uma fungao teste. 

Convem observar que neste trabalho a parte temporal de (5.281 e discretizada 
utilizando os aproximantes de Fade, Rn e R 22 , e a parte espacial via formulagoes 
dos metodos de elementos finitos: MEFMQ, MEFG e SUPG. 


6. Resultados Numericos 

Todos os resultados apresentados na sequencia resultam das formula^oes semi- 
discretas abordadas. Apresentamos analises dos erros numericos a partir das normas 
L 2 e Loo, comparando as soluQoes numericas com as solugoes anah'ticas dos exemplos 
avaliados. 


6.1. Exemplo 1 

A equagao ID de Burgers definida em ( |2.1[ ) com condigao inicial 

2 e 7 r sin( 7 ra;) 


u{x, 0) = 


fc > 1, 


k + cos(7rx) ’ 

e condigoes de fronteira u{0,t) = 0 = u{l,t) tern como solugao analitica [TB] 

2e7r exp(7r^et) sin(7ra:) 


i{x, t) = 


k + exp(— TT^et) cos(7ra:)' 


com fc > 1. 


(6.29) 


(6.30) 


Gonsiderando o dominio 0<a;<lefc = 2, apresentamos nas Figuras la-lf 
os resultados das formulagoes semi-discretas comparados com o resultado anah'tico, 
para os tempos t = 0.5 e t = 1, com At = 0.5, Re = 10® e uma malha de 50 
elementos lineares, o que equivale a h = 0.02. 

Observamos que o aproximante de Fade Ru adicionado as formulagoes MEFMQ, 
MEFG e SUPG, apresentou oscilagoes mais evidentes no tempo final t = 1 con- 
forme Figuras la-lc, mas ressaltamos a ocorrencia de menor intensidade para o caso 
MEFMQ+i?ii, Figura la. Ao incluirmos o aproximante de Fade R 22 nas formu¬ 
lagoes MEFMQ, MEFG e SUPG, verificamos que as oscilagoes foram amenizadas, 
Figuras Id-lf. 

Para uma avaligao mais precisa das formulagoes semi-discretas, calculamos os 
logaritmos dos erros considerando malhas com h = 1/50 e h = 1/1000 para os 
tempos t = 0.5 e t = 1 com At = 0.5 e Re = 10®, tais resultados encontram-se na 
Tabela 2. 

A partir dos resultados apresentados na Tabela 2, confirmamos que o aproxi¬ 
mante de Fade R 22 em conjunto com as formulagoes MEFMQ, MEFG e SUPG 
aumentou a regiao de convergencia das solugoes numericas e apresentou maior pre- 
cisao quando comparado as solugoes obtidas por meio de aproximantes de Fade 
i?ii, considerando h = 1/50. Melhorando o refinamento na diregao espacial, isto, e, 
tomando h = 1/1000, observamos que a formulagao MEFMQ+i?ii mostrou-se mais 
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Figura 1: a-f) Solugao anali'tica e soluQoes numericas das formulagoes semi-discretas. 


Tabela 2: Calculos dos logaritmos dos erros entre as formulagoes semi-discretas 
utilizando as normas L 2 e Loo- 




h = 

1/50 

h = 

1/1000 

Formulagoes E(logiQ) 

t = 0.5 

t = 1 

t = 0.5 

t = 1 

MEFMQ+Lii 


-1.4338 

-1.0964 

-1.4709 

-1.1762 


L||oo 

-1.0570 

-6.6807x10-^ 

-1.2759 

-9.7797x10-1 

MEFMQ+L 22 

\Eh 

-1.4363 

-1.0998 

-1.4649 

-1.1698 


E\U 

-1.0562 

-6.4226x10"^ 

-1.2566 

-9.5395x10-1 

MEFG+i?ii 

\Eh 

-1.4117 

-8.6034x10"^ 

-1.4614 

-1.1665 


E\\oo 

-1.0112 

-5.1914x10-1 

-1.1212 

-8.1722x10-1 

MEFG+i?22 

\Eh 

-1.4289 

-1.0365 

-1.4612 

-1.1664 


L||oo 

-1.0261 

-5.8131x10-1 

-1.1178 

-8.1754x10-1 

SUPG+i?ii 

|£^I|2 

-1.4130 

-8.7854x10-1 

-1.4614 

-1.1665 


L||oo 

-1.0128 

-5.2995x10-1 

-1.1212 

-8.1722x10-1 

SUPG+i?22 

|£^I|2 

-1.4286 

-1.0351 

-1.4615 

-1.1665 


L||oo 

-1.0260 

-5.8146x10-1 

-1.1514 

-8.3593x10-1 


precisa entre todas as formulagoes, nao havendo entao a necessidade de um metodo 
oneroso para a discretizagao no tempo. 
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6.2. Exemplo 2 

Seja um problema de propagagao uniforme de cheque j22| . para a equagao ID de 
Burgers (2.11, com condigao inicial dada por 


u{x, 0) = 
cuja solugao anah'tica e dada por 
u{x, t) = Ui + 


Ml = 0.5, se X < 0 
U 2 = 1.5, se a: > 0, 


M2 — Ml 

1 + exp [-Re^i^ix - 


(6.31) 


(6.32) 


Considerando o domi'nio —0.5 < x < 0.5 e as condigoes de fronteira satisfazendo 
m(— 0.5,0) = 0.5 e m(0.5,0) = 1.5, apresentamos nas Figuras 2a-2f os resultados das 
formulagoes semi-discretas comparados com o resultado analitico, para os tempos 
t = 0.05 et = O.l, com At = 3.3 x 10“"^, Re = 10"^ e uma malha com 3000 elementos 
lineares. 


a) <>) c) 









• MEFMQ, flu, 1 = 0.05 

2 

^ MEFGiflu. 

II 

o 

© 


SUPG 1 /?ii. 

= 0.05 


0 . . MEFMQ 1 Ri 

. t=Q.l 


® MEFG-/?[,. 

= 0.1 


0 SUPG • /?!,. 

f =0.1 

rs 

-Analitica 


-Analftica 


-Analitica 

1.6 


1.6 




1.6 


1.4 




1.4 




1.4 




■S 




1.2 




1.2 




1 




1 












0.8 




0.8 




0.6 




0.6 




0.6 




0.4 




0.4 




0.4 




0.2 


0.2 


0.2 



-0.5 0 0.5 ^.5 0 0.5 .0.5 o 0.5 





Figura 2: a-f) Solugoes numericas das formulagoes semi-discretas para h = 1/3000, 
At = 3.3 X 10“'*, Re = 10^ e para os tempos t = 0.05 e t = 0.1. 

Podemos observar nas Figuras 2d-2f que ao adicionarmos o aproximante de Fade 
i ?22 nas formulagoes MEFMQ, MEFG e SUPG, obtivemos um aumento na regiao de 
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convergencia das solugoes numericas, assim como uma maior precisao das solugoes 
quando comparadas as obtidas por meio do aproximante de Fade Rn, vide Figuras 
2a-2c. 


Tabela 3: Calculos dos logaritmos dos erros entre as formulagoes semi-discretas 
utilizando as normas L 2 e Loo 

h = 1/3000 /i= 1/4000 


Formulagoes E(logig) t = 0.05 t = 0.1 t = 0.05 t = 0.1 


MEFMQ+Lii 

\E\\2 -1.9468 -1.8369 -1.8382 -1.7146 

LI 00 -6.5205x10-1 -5.8740x10-1 -6.4089x10-1 -5.9775x10-1 

MEFMQ+L 22 

\E\\2 -2.1093 -1.9833 -1.8465 -1.7192 

LI 00 -8.6062x10-1 -7.7387x10-1 -6.8386x10-1 -6.2629x10-1 

MEFG+Lii 

|L|2 -1.9063 -1.8109 -1.8397 -1.7156 

L| 00 -6.2961x10-1 -5.7230x10-1 -6.3991x10-1 -5.9758x10-1 

MEFG+L 22 

|L|2 -2.2371 -2.0699 -1.8497 -1.7206 

L| 00 -9.1682x10-1 -8.4416x10-1 -6.8569x10-1 -6.2722x10-1 

SUPG+i?ii 

|L|2 -1.9063 -1.8109 -1.8397 -1.7156 

L| 00 -6.2961x10-1 -5.7230x10-1 -6.3991x10-1 -5.9758x10-1 

SUPG+L 22 

|L|2 -2.2473 -2.0729 -1.8496 -1.7206 

L| 00 -9.4179x10-1 -8.5595x10-1 -6.8579x10-1 -6.2725x10-1 


Porem, confrontando MEFMQ+L 22 com as formulagoes MEFG+L 22 e 
SUPG+i ?22 observamos ainda a ocorrencia de pequenas oscilagoes na vizinhanga 
do choque, Figura 2d, quando h = 1/3000. Estas oscilagoes persistem quando 
utilizamos h = 1/4000, como pode ser confirmado observando a Tabela 3, onde 
apresentamos os calculos dos logaritmos dos erros entre as formulagoes estudadas. 
Uma alterantiva para amenizar estas oscilagoes seria utilizar um estabilizador no 
MEFMQ. 

Por fim, verificamos que ao utilizarmos At da ordem de 10“^ as formulagoes 
MEFG+i ?22 e SUPG+i ?22 possibilitaram resultados mais apurados quando compa¬ 
radas as demais formulagoes semi-discretas, independente do refinamento do espago. 


7. Conclusao 

O trabalho forneceu uma comparagao entre as formulagoes semi-discretas para re¬ 
solver a equagao ID de Burgers, para diferentes condigoes iniciais e de fronteira. 
Observamos que o aproximante de Fade R 22 adicionado as formulagoes MEFMQ, 
MEFG e SUPG apresentou solugoes aproximadas com convergencia mais rapida e 
maior precisao em comparagao ao i?ii. Tambem, o metodo R 22 amenizou conside- 
ravelmente as oscilagoes quando usado nas formulagoes MEFG e SUPG. 

Verificamos no Exemplo 1 que melhorando o refinamento na diregao espacial a 
formulagao MEFMQ+i?ii mostrou-se mais precisa entre as demais formulagoes, nao 
havendo a necessidade de um metodo oneroso para a discretizagao no tempo. En- 
quanto que no Exemplo 2 as formulagoes MEFG+i ?22 e SUPG+L 22 possibilitaram 
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resultados mais apurados quando comparadas as outras formulagoes semi-discretas, 
independente do refinamento do espago. Salientamos entao que, para os exemplos 
apresentados, predominou a escolha da formulagao para a discretizagao temporal e 
espacial. 

Notamos ainda que os erros obtidos nas formulagoes semi-discretas, ao refinar- 
mos a malha, diminuiram consideravelmente, apresentando concordancia entre as 
solugoes numericas com a solugao exata independente da formulagao utilizada. 

Abstract. In this work we compare semi-discrete formulations to obtain nu¬ 
merical solutions for the ID Burgers equation. The formulations consist in the 
discretization of the time-domain via multi-stage methods of second and fourth 
order: i?ii and R 22 Bade approximants, and of the spatial-domain via finite ele¬ 
ment methods: least-squares (MEFMQ), Galerkin (MEFG) and Streamline-Upwind 
Petrov-Galerkin (SUPG). Knowing the analytical solutions of the ID Burgues equa¬ 
tion, for different initial and boundary conditions, analyzes were performed for nu¬ 
merical errors from L 2 and Loo norm. We found that the R 22 Fade approximants, 
added to the MEFMQ, MEFG, and SUPG formulations, increased the region of 
convergence of the numerical solutions, and showed greater accuracy when compa¬ 
red to the solutions obtained by the Ri\ Fade approximants. We note that the R 22 
Fade approximants softened the oscillations of the numerical solutions associated 
to the MEFG and SUPG formulations. 
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